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We investigate the effects of an external magnetic helicity production on the evolution of the 
cosmic axion held. It is shown that a helicity larger than {few x 10“^®G)^Mpc, if produced at 
temperatures above a few GeV, is in contradiction with the existence of the axion, since it would 
produce too much of an axion relic abundance. 
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Recently, the topic of cosmic magnetic helicity genera¬ 
tion and phenomenolo^ has been widely discussed (see, 
e-g-, An in depth study of this 

very peculiar quantity could help in understanding the 
nature of the cosmic magnetic field itself and its genera¬ 
tion. Magnetic helicity is defined as 

HB{t) = d^xA-V X A, (1) 

where = (A°,A) is the electromagnetic field. In a 
flat universe described by a Robertson-Walker metric, 
where R{t) is the expansion parameter 
normalized so that at the present time to, R{to) = 1, the 
electric and magnetic fields are defined as E = —R~^A 
and B = x A, where a dot indicates the derivative 

with respect to the cosmic time t. Hence, the helicity can 
be written as Hsit) = {1/V)R^ jy d^x A ■ B. It should 
be noted that in the literature the definition oi Hb is 
usually given without the factor i?^. In any case, the 
two definitions coincide at the present time. However we 
note that, since the early universe is a very good conduc¬ 
tor, magnetic fields are frozen into the plasma and then 
B scales in time as i? oc R~^ Q. Therefore, magnetic 
helicity, as defined in Eq. o, remains constant (after its 
generation). 

The magnetic helicity is related to the topological 
properties of the magnetic field, it is a CP—odd pseudo¬ 
scalar quantity and, if different from zero, would reveal 
a macroscopic P and CP violation in the universe. 

Let us introduce the magnetic energy and magnetic he¬ 
licity spectra £B{k,t) = (27r/R)A:^ B(k) • B(—k) and 
Ti.B{k,t) = (47r/R)P^fc^A(k)-B(—k), respectively, where 
B(k) and A(k) are the magnetic field and the vec¬ 
tor potential in Fourier space and k = |k|. In 

terms of the spectra, the magnetic energy, EB{t) = 
{\/2V)jyd^X^^ (x, t) , and the magnetic helicity are 
Eb = jdk£B, and Hb = JdkHB- From the above 
definitions it follows that any magnetic held conhgura- 
tion satishes the inequality \HB{k,t)\ < 2R‘^k~^£B{k,t) 
1^. A straightforward integration in k leads to the 
so-called “realizability” condition, \Hb\ < = 

2,R?£,bEb, which hxes the maximal helicity associated 
with a given magnetic held conhguration. Here, f,B{t) = 


2ttR f dk k~^£B{k, t)/EB is the magnetic held correlation 
length. Therefore, bounds on the magnetic held can be 
directly translated into bounds on magnetic helicity. Be¬ 
fore discussing this point more deeply, it is useful to intro¬ 
duce the following parametrization for the helicity pro¬ 
duced at the temperature T*: 

Pb(A) = ±3.3 X 10-1® r 9^ G'Mpc, (2) 

where 5 * and g^,s count the total numbers of effectively 
massless degrees of freedom referring to the energy and 
entropy density of the universe, respectively. (From now 
on, g* and ( 7*5 will be considered equal.) The choice r = 1 
gives the maximal helicity associated with a magnetic 
held correlated on the Hubble scale, ^b = H~^ {H is the 
Hubble parameter), and whose energy density equals the 
total radiation energy density at T*. Thus, the parameter 
r < 1 measures the fraction of helicity with respect to the 
maximum possible at a certain time. 

It is known from the Big Bang Nucleosynthesis (BBN) 
argument that a magnetic held B at the temperature 
T ~ O.lMeV, correlated on the Hubble scale at that time, 
cannot be more intense than IQHG 0. This yields the 
upper limit r < Tben — 82 ( 7 *(rH.)i/^ (T*/GeV). On the 
other hand, the analysis of the observed anisotropy of 
the Gosmic Microwave Background (CMB) radiation re¬ 
quires the bound 10-®G on the scales of IMpc 0 for 
the magnetic held (today), which leads to r < Tcmb — 
0.2 ( 7 h.(T*)^/^ (P,/GeV). Observe that for magnetic helic¬ 
ity produced before IGeV neither of these limits give any 
interesting constraints, since r is always expected to be 
less than 1 . 

In this paper we consider the effects of magnetic he¬ 
licity generation on the cosmic axion held, showing how 
such generation could considerably amplify the axion’s 
expected relic abundance £la- In particular, we will show 
that the request £la < ^matter sets, in most cases, a 
bound on r of many orders of magnitude smaller that 
the above Tebn and Tcmb- 

The axion held a emerges in the so called Pecce-Quinn 
(PQ) mechanism 0 for the solution of the strong GP 
problem (for a review see, e.g., M)- In most models, 
it is identihed with the phase 0 of a complex scalar held 
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a = faQ, where fa is a phenomenological scale usually re¬ 
ferred to as the PQ- or axion-constant, and presently con¬ 
strained in the very narrow region 10 ® ^ fa lO^^GeV 
by astrophysical and cosmological considerations 0 . 1 
For temperatures above the PQ-scale the total La- 
grangian is required to be invariant under a constant shift 
of the phase 0 (PQ-symmetry). However, for lower tem¬ 
peratures, the symmetry is (spontaneously) broken and 
the angle 0 is fixed on a precise value. We indicate this 
initial value with 0^. This phase evolves following the 
equation ^ 


•• • (JV 

0 + 3iJ0 + = 0, (3) 

where F(0,T) = m®(T)(l — cos0) is the instanton in¬ 
duced potential. The temperature-dependent axion mass 
is 0, ^ ma{T) ~ O.lm (A/T)^'^ for T ^ A and 
ma{T) = m for T ^ A with m ~ 6.2 x 10 “®eV//i 2 
[/i 2 = /a/(10^®GeV)], and A^200MeV is the QCD scale. 

As the curvature (mass) term in Eq. becomes domi¬ 
nant over the friction (Hubble) term, 0 begins to oscillate 
with the frequency ma{T). This happens at about the 
temperature T = Ti defined by the equation maiTi) = 
?>H{Ti). Approximately Ti ~ 0.9/]^®'^’^^A25o® GeV with 
A 200 = A/(200MeV). 

During this period of coherent oscillations the num¬ 
ber of axions in a comoving volume remains constant, 
UaB? = const, where na(T) is the axion number density. 
Thus, the axion relic abundance today is 

~ 0.2 A 2 -o°o®" [0? + (0i/3Afi)" ] (4) 

where 0i, 0i and Hi are respectively 0(T'i), ©(Ti) and 
H{Ti). If the axion field is not interacting in the region 
T > Ti, and Ti fa, then the axion kinetic energy 
at T = Ti is reliably negligible, 0 i oc (Ti/Za)^ 'C 1 , 
and consequently 0i ~ 0^. Thus, Eq. 0 reduces to 
= ^ao — 0-2 AAo®^ 0 |/i 2 ^^^- Assuming the natural 
value Qi ~ 1, this gives ilao — 0.3 (the expected value 
for the dark matter abundance) for /12 ~ 1. Much larger 
values of /12 would cause too much axion production and 
are therefore excluded. 

However, if the field evolution took place in the pres¬ 
ence of an external magnetic field, the previous analysis 
needs to be modified. It is easy to show the importance 
of the magnetic helicity in this case. In fact, in the pres¬ 
ence of a magnetic field, the right hand side of Eq. © is 
modified as (a/27r/^)E • B. Assuming a homogeneous 
axion field, and observing that (I/F) fyd^x E B = 


-HbI{2R^), Eq. © becomes 


dV 

0 + 3i70 + — 
a0 


a Hb 
4^/2 i ?3 ■ 


(5) 


The slightness of the effects of a non-helical magnetic field 
on the cosmic evolution of the axion field discussed in 
literature (see, e.g., ^ 3 ) can easily be understood: even 
though the axion itself produces helicity in the presence 
of an external magnetic field [3,13 and therefore its 
evolution is modified, the effect is negligible since the 
helicity produced by the axion itself is small. However, if 
another mechanism were responsible for the production 
of a significant amount of helicity, the consequences on 
the axion relic abundance today could be enormous. 

Let us assume, first, that a large amount of helicity is 
produced at the temperature T, ^ Ti and that the time 
scale of its production is short compared to the time scale 
ti. ® Moreover, taking into account that after its gener¬ 
ation the helicity can be considered a (quasi)-conserved 
quantity (see, e.g., jl^l), we can safely approximate it as 
a step function, Lls(t) = iL_B(G) S{t — t*). 

When the potential energy is small compared to the 
kinetic energy and/or the friction, Eq. (0 can be solved 
giving 0(r) = aHB{T)/{ATTffE?), that is, the kinetic 
energy is a decreasing function of time. Let us indicate 
with Teq the temperature at which the kinetic energy 
equals the potential: (0^/2)|Teq = (^(0))|Teqj where for 
the sake of simplicity we are considering the root mean 
square value of the potential (E(0)) = m®(T)(I —1/-\/2). 
It is clear that only if Teq > Ti will the oscillations start 
at the temperature Ti , where the friction becomes small. 
In this case, the standard analysis applies and we find 
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1 + {r/raf 


if ?’ < ra, ( 6 ) 


where r < ~ I.5x I 0 -“( 7 y"(TAAAo®Vi¥^®(ir*/GeV) 

corresponds to Teq > A. However, if Tgq < A the oscil¬ 
lations cannot start at T = A) the kinetic energy being 
still too large, but will be delayed and start at T = Aq. 
The resulting axion relic abundance is 

Da ~ 0.2AV^®A^o*cf (r/ra), if r > r^. (7) 

In Fig. Q] we show the cosmologically allowed region, 
Da < Dniatter — 0.3, for the PQ-constant A and the 
r parameter. Observe that, since we are now consider¬ 
ing helicity to be generated before the axion oscillations, 
A > A ~ IGeV, the BBN and GMB arguments do not 
set any consistent limit on the magnetic helicity. On the 
other hand, the limit on r from the axion abundance is 


^ For the so-called hadronic axion models, whose interaction with 
electrons is suppressed, a small region around fa ~ 10®GeV is 
also allowed For other possibilities of re cons idering the 

limits on the PQ-constant see, for example, Ref. 0. 


^ This is not a very strong assumption. For example, in the a- 
dynamo mechanism this time scale is expected to be At* ^ "C 

ft 0 . 
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FIG. 1: Allowed region for the magnetic helicity and the PQ- 
constant in the case of helicity production before the axion 
coherent oscillations. The area under the curve represents the 
parameter region Qa < ^matter ^0.3 m. 


extremely significant: the curve in Fig. ^ is well approx¬ 
imated by r = Ta (for /12 < 1) and this implies [see 
Eq. 1(2)] that a magnetic helicity generated at the tem¬ 
perature r* cannot be larger than 

Hb^5x 10-30 A-o- 03/22175 Q^Mpc. (8) 

In the most conservative case, /12 ~ 1, we exclude any 
Hb > (few X 10-30G)^Mpc. 

We note that a full numerical treatment of the axion 
evolution equation would reveal the existence of a small 
region with fa above the cosmological bound. We omit 
the discussion of this case to avoid useless complications. 
Whatever the case, the region emerges as a fine tuned 
compromise between helicity and the 0-angle, and is 
therefore very small. 

We will now briefly discuss the phenomenon of helicity 
production during the axion oscillations and show that 
the effects on the axion relic abundance are, in this case, 
less evident. The peculiarity of this second case is the ap¬ 
pearance of an additional energy scale, the axion mass, 
which was neglected in the above discussion. In the radi¬ 
ation era Eq. 0 can be re-stated as /4:+ma)ip = 

F{t), where (p = 0i?3/2^ F{t) = -(a/47r/^)(ijB/i?3/2)^ 
and we assumed |0| < I (see below) so that L(0) — 
rriaQ^f^. For temperatures T <C A the axion mass is con¬ 
stant, mair) = m, and much larger then the Hubble pa¬ 
rameter, H/m ~ 10-®A2 oo/i2 5 *('E)(T/A)^. Hence, fac¬ 
tors of order H/rua can be neglected in the previous equa¬ 
tion, which therefore describes a forced harmonic oscilla¬ 
tor with frequency m and forcing term F{t). Its solution 
is (^ = (/jo + (1/to) J^_duF{u) sinm(t —m), where ipo is the 
solution for the free harmonic oscillator. With a suitable 
choice of initial conditions we can write (po = Asin(mt), 
A being related to the actual axion energy density in the 
case of null magnetic helicity, Pa{To) = ^famfA^. 


Let us suppose that the helicity is produced during the fi¬ 
nite interval of time {Ia-, ts)- Since before Ia and after Ib 
the evolution of p is that of a free oscillator, and since for 
f ^ ?'cMB it results |0| < 1 if T* < IMeV, we can apply 
the standard procedure to find the axion relic abundance 
today. This reads fla = + 2H^Q^Rez -I- jzp, where 

■z = {falV‘^Pc) Jt^dtF{t)e^‘^*, Pc is the present criti¬ 
cal density, and Clao is the expected axion relic abun¬ 
dance (today) in the absence of magnetic helicity [see 
the discussion below Eq. 0 ]. The time profile of the 
magnetic helicity in the expanding universe is conve¬ 
niently parameterized by the following “rapidity” func¬ 
tion, /(<) = (i?*/ii)3/2 G Hb/F[B{ t»), where L is a refer¬ 
ence time between Ia and Ib- Therefore, the forcing term 
can be written as F{t) = — (a/ 47 r/^)[iJB(T*)/ii*^^]/(t) 
and, using Eq. 0 , z becomes 

z - ±3 X /-V (r*/A)i/2 /(m), (9) 

where /(w) = dtf{t) is the Fourier transform 
of f{t). The axion relic abundance can be exactly cal¬ 
culated once the rapidity function f[t) is known. In 
the case of slow (At ^ m~^) helicity production it re¬ 
sults f{m) ~ 0 and, consequently, fla — ^ao- On the 
other hand, if helicity production is fast (At <C m~^) 
we have |/(m)| ~ Re[/(m)] ~ 1. To give a simple esti¬ 
mate of ila we observe that for /12 < 1 it results that 
flao ^ fla) and therefore fla ^ \z\'^. This yields the up¬ 
per bound r < IQ-^A^qq /12 (A/T*)^/^ which translates 
into a bound on the present magnetic helicity 

Lis < 10-2"A3/2 5 *(r,)-i/Vi 2 (^) ' G^Mpc. (10) 
This limit is always stronger than the one from the BBN 

2 /3 

and, for T* > few fy^ keV, it also results that r < Tcmb 
( see Fig. 0 . In this figure, we have included the region 
above the cosmological limit fa > lO^^GeV (see the con¬ 
tinuous line), since the analysis is straightforward in this 
case. Simply note that the minus sign in Eq. © corre¬ 
sponds to a reduction of the axion relic abundance. 

To conclude, we have pointed out that a strong con¬ 
nection between magnetic helicity and axion cosmology 
exists. In particular, we have studied the influence of pri¬ 
mordial magnetic helicity generation on the present axion 
relic abundance. The result is that, allowing for possi¬ 
ble helicity generation in the early universe, the knowl¬ 
edge of the PQ-constant would not be sufficient to de¬ 
termine the axion relic abundance today. Indeed, sup¬ 
pose axions were detected with the PQ-constant strictly 
under lO^^GeV. This would not exclude the possibil¬ 
ity that they represent the dark matter component of 
our universe, if an appropriate amount of helicity has 
been generated (see Fig. 1). We focused our analy¬ 
sis on two different cases corresponding to helicity pro¬ 
duction before or during the axion coherent oscillations. 
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FIG. 2: Allowed region for the magnetic helicity and the 
PQ-constant in the case of fast helicity production during 
the axion coherent oscillations. The areas enclosed by the 
dashed and continuous curves represent the parameter region 
< fimatter — 0.3 for the plus and minus sign in Eq. 
respectively. The areas under the horizontal lines represent 
the region allowed by the CMB analysis for three different 
temperatures. 


Let us consider the former, and more interesting, case 
first. We look at two different possibilities within this 
case: i) first, if axions were detected as cold dark mat¬ 
ter particles this would severely constrain the allowed 
amount of magnetic helicity in our universe, according 
to Eq. 0. On the other hand, ii) if a helicity greater 
than {few x 10“^®G)^Mpc were revealed at the present 
time this would be in contradiction with the existence 
of the invisible axion itself, unless it is assumed that he¬ 
licity production took place at a temperature above the 
PQ-scale. It should be noted that the bounds found are 
considerably stronger than the ones deduced from the 
BBN or CMB analysis. 

Finally, if helicity production took place during the 
primordial axion oscillations, the constraints on the he¬ 
licity depend on the time scale of its production, and are 
indeed relevant only if this is much smaller than the typ¬ 
ical axion oscillation time 1 /toq. Also in this case, the 
bounds found, Fig. Inland Eq. cni), are relevant compared 
to the ones from the BBN and CMB. 

We consider these results of interest in view of next 
generation experiments devoted to the axion search and 
future experiments on CMB, which are expected to be 
sensitive enou^ for the detection of a possible helical 
magnetic field^. Ei]|. 
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